ftr®G^D@[r@ Don 

PLANT SCIENCE 



ORIGINAL RESEARCH ARTICLE 

published: 15 March 2013 
doi: 10.3389/fpls.2013. 00050 




Grass cell wall feruloylation: distribution of bound ferulate 
and candidate gene expression in Bra chy podium 
distachyon 

Hugo B. C. Molinari 1 ', Till K. Pellny, Jackie Freeman, Peter R. Shewry and Rowan A. C. Mitchell* 

Plant Biology and Crop Science, Rothamsted Research, Harpenden, Hertfordshire, UK 



Edited by: 

Maurice Bosch, Institute of Biological, 
Environmental and Rural Sciences, 
Aberystwyth University, UK 

Reviewed by: 

Ana Winters, Aberystwyth University, 
UK 

Carsten Rautengarten, Joint 
BioEnergy Institute, Lawrence 
Berkeley National Laboratory, USA 
Jaime Barros-Rios, Umea University, 
Sweden 

'Correspondence: 

Rowan A. C. Mitchell, Plant Biology 
and Crop Science, Rothamsted 
Research, Harpenden, Hertfordshire 
AL5 2JQ, UK. 

e-mail: rowan. mitchell@rothamsted. 
ac.uk 

1 Present address: 

Hugo B. C. Molinari, Laboratory of 
Genetics and Biotechnology, Embrapa 
Agroenergia, Brasilia, Federal District 
70770-901, Brazil. 



The cell walls of grasses such as wheat, maize, rice, and sugar cane, contain large amounts 
of ferulate that is ester-linked to the cell wall polysaccharide glucuronoarabinoxylan (GAX). 
This ferulate is considered to limit the digestibility of polysaccharide in grass biomass as it 
forms covalent linkages between polysaccharide and lignin components. Candidate genes 
within a grass-specific clade of the BAHD acyl-coA transferase superfamily have been 
identified as being responsible for the ester linkage of ferulate to GAX. Manipulation of 
these BAHD genes may therefore be a biotechnological target for increasing efficiency 
of conversion of grass biomass into biofuel. Here, we describe the expression of these 
candidate genes and amounts of bound ferulate from various tissues and developmental 
stages of the model grass Brachypodium distachyon. BAHD candidate transcripts and 
significant amounts of bound ferulate were present in every tissue and developmental 
stage. We hypothesize that BAHD candidate genes similar to the recently described Oryza 
sativa p-coumarate monolignol transferase (OsPMT) gene (PMT sub-clade) are principally 
responsible for the bound para-coumaric acid (pCA), and that other BAHD candidates (non- 
PMT sub-clade) are responsible for bound ferulic acid (FA). There were some similarities 
with between the ratio of expression non-PMT/PMT genes and the ratio of bound FA/pCA 
between tissue types, compatible with this hypothesis. However, much further work to 
modify BAHD genes in grasses and to characterize the heterologously expressed proteins 
is required to demonstrate their function. 

Keywords: glucuronoarabinoxylan, hydroxycinnamic acid, BAHD gene family, PF02458 domain, bound phenolic 



INTRODUCTION 

Many of the most abundant potential lignocellulosic feedstocks 
are from grasses, whether by-products from food crops such as 
wheat straw, rice straw, maize stover, and sugar cane residues or 
from specialized bioenergy crops such as Miscanthus and switch 
grass. Such feedstocks differ from those from dicotyledonous crops 
in the occurrence of frans-ferulic acid (FA), a hydroxycinnamic 
acid which is ester-linked to the hemicellulosic component, glu- 
curonoarabinoxylan (GAX; Carpita, 1996; Scheller and Ulvskov, 
2010). This linkage occurs only in commelinid monocots and is 
particularly abundant in grasses (Harris and Trethewey, 2010). It 
is considered to affect digestibility (Grabber, 2005; Dhugga, 2007), 
as FA can oxidatively cross-link to form ether bonds or C-C bonds, 
linking chains of GAX (Ishii, 1991), or of GAX to lignin (Ralph 
etal., 1995). Thus it covalently links the cell wall polysaccharide 
which can be used as a substrate for biofuel to the inhibitory lignin 
component. Whilst the barrier to degradation presented by FA is 
mainly structural, FA also has antimicrobial activity, inhibiting 
sugar fermentation by yeast (Baranowski et al, 1980; Akin, 2008). 
Grasses with decreased feruloylation of GAX in their cell walls 
may therefore be more readily converted into biofuel. In partic- 
ular, they will have a lower energy requirement for separation of 
lignin from polysaccharide, an essential step in most pipelines for 
biomass processing. 



The genes and enzymes responsible for feruloylation of GAX 
are still unknown. However, using a bioinformatics analysis of 
publicly available expressed sequence tags (ESTs), we identified 
an orthologous group within the BAHD acyl-coA superfamily 
as potential candidates based on their differential expression in 
cereals compared with dicots (Mitchell etal., 2007). The BAHD 
family (also referred to as the PF02458 family, as all family mem- 
bers contain the PFAM domain PF02458 which is specific to 
the family) is named after the first four members to be bio- 
chemically characterized (benzylalcohol acetyltransferase, BEAT; 
anthocyanin hydroxycinnamoyl transferase, AHCT; anthranilate 
hydroxycinnamoyl/benzoyl transferase, HCBT; deactylvindoline 
acetyltransferase, DAT; D'Auria, 2006). The predicted role in GAX 
feruloylation is compatible with the several hydroxycinnamoyl 
transferase activities creating ester bonds known in the family, 
with feruloyl-CoA acting as the donor. We also identified GT61 
genes as likely to be involved in GAX biosynthesis (Mitchell et al, 
2007), and have recently demonstrated (in collaboration with 
the Dupree lab, University of Cambridge, UK) that these genes 
mediate the addition of 3-linked arabinofuranose (Ara/) to xylan 
(Anders et al., 2012). GT61 genes tend to be co-expressed with the 
BAHD candidates in cereals (Mitchell etal., 2007; Shewry etal, 
201 1), consistent with a role in feruloylation of the 3-linked Ara/ 
residues added by GT61 proteins. Our predicted role of BAHD 
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genes in feruloylation was tested by another group using RNAi 
suppression in rice, showing a modest (—19%) but significant 
decrease in cell wall bound ferulate in the stems of the transgenic 
plants (Piston etal, 2010). However, this effect could be an indi- 
rect consequence of gene suppression, so the role of these BAHD 
genes still requires elucidation. 

One question is the intracellular localization of the proteins 
encoded by the BAHD genes: this is predicted to be cytosolic based 
on the protein sequences and this is the localization for all known 
members of the family (D'Auria, 2006), whereas GAX is synthe- 
sized in the Golgi (Buanafina, 2009). An activity from rice seedlings 
capable of feruloylation of Ara/on a synthetic molecule was found 
in the soluble fraction rather than the membrane fraction, con- 
sistent with a cytosolic enzyme (Yoshida-Shimokawa et al., 2001). 
A possible explanation is that BAHD proteins are responsible for 
feruloylation of a cytosolic precursor, such as UDP-Ara/ which is 
the substrate for xylan arabinosylation. Whereas the great major- 
ity of UDP-arabinopyranose is synthesized in the Golgi lumen, 
this is converted by a mutase to the UDP-Ara/ form in the cytosol 
(Konishi et al., 2007; Rautengarten et al., 20 1 1 ) so this is a possible 
substrate for a cytosolic feruloylation reaction inplanta. The prod- 
uct from this would then be transported back into the Golgi, where 
feruloylated Ara/ would be transferred onto GAX, possibly by a 
GT61 enzyme. The feruloylation of arabinoxylan has been shown 
to occur within the protoplast (Myton and Fry, 1994), although 
there is also evidence of an alternative pathway which allows fer- 
uloylation in the wall when inhibitors are used to suppress the 
secretory pathway (Mastrangelo et al., 2009). 

A key recent breakthrough was the demonstration that at least 
one of these BAHD genes in rice has a different role from that 
predicted; in that the encoded protein catalyzed the addition 
not of FA to a sugar, but of the closely related para-coumaric 
acid (pCA) to monolignols. This gene therefore encodes a pCA 
monolignol transferase and has been named OsPMT ( Oryza sativa 
p-coumarate monolignol transferase; Withers etal, 2012). Inter- 
estingly, whereas other BAHD proteins have been shown to use 
both FA-CoA and pCA-CoA as donor molecules (Luo et al., 2007), 
OsPMT was found to have almost no activity with FA-CoA (With- 
ers etal, 2012). Whereas it is theoretically possible that OsPMT 
could also add FA or pCA to Araf, activities with such different 
acceptors have not previously been found for BAHD proteins. 
OsPMT activity is presumably responsible for forming ester link- 
ages between pCA and lignin, a linkage which is also much more 
common in grasses than in other groups of plants. It is there- 
fore possible that the candidate group which includes OsPMT is 
responsible only for this activity, rather than also for GAX fer- 
uloylation. However, to date, only one of the genes has been 
shown to have this activity and others are expressed (albeit at 
low levels) in wheat starchy endosperm which does not have 
lignin but does have feruloylated AX (Pellny etal, 2012), sug- 
gesting that other genes within the group are responsible for 
xylan feruloylation. During preparation of this article, strong 
support for this theory was published where another BAHD 
gene within this group was overexpressed in rice, resulting in 
increased bound pCA with good evidence that this was ester- 
linked to GAX, not lignin (Bartley etal, 2013). Whereas this 
gene seems specific for pCA, it seems very likely that some of the 



other similar genes in the clade are responsible for FA ester-linked 
to GAX. 

Brachypodium distachyon is a model grass species which has a 
small and fully sequenced genome, a short life cycle, and small size, 
making it ideal for the study of grass cell walls (Vogel et al., 2010). 
We have therefore determined the contents of bound FA and pCA 
(predominantly linked to GAX and lignin, respectively) in various 
tissues and at different stages throughout the life cycle, and related 
these to the expression of BAHD genes identified as candidates for 
the addition of these phenolic acids to cell wall polymers (GAX 
and lignin). 

RESULTS AND DISCUSSION 

Figure 1 shows a phylogenetic tree of the clades of the BAHD gene 
superfamily which have been proposed to contain candidate genes 
for the feruloylation of GAX in grass cell walls (Mitchell etal., 
2007), divided into clades A and B. Clade B genes are not highly 
expressed or co-expressed with xylan pathway genes, and do not 
have clear orthologs between rice and Brachypodium. The more 
highly expressed clade A genes were therefore identified as the 
stronger candidates (Mitchell etal, 2007) for AX feruloylation. 
Clade A also includes the gene subsequently shown to encode a 
monolignol pCA transferase in rice, OsPMT (Withers et al., 20 12). 
Two wheat genes, orthologous to candidates 1 and 3, are expressed 
in wheat starchy endosperm, which does not contain lignin in the 
cell walls (Pellny etal., 2012). An RNAi construct which simul- 
taneously suppressed the expression of rice candidates 1, 4, 8, 
and 10 (in our nomenclature) resulted in a decrease in cell wall 
FA of 19% in transgenic rice (Piston etal., 2010). A possible 
interpretation is therefore that the monolignol pCA transferase 
activity is restricted to the sub-clade containing OsPMT which 
includes candidates 7 and 9 and the only Arabidopsis gene, and the 
other candidates in clade A may be involved in GAX feruloylation. 
This is supported by the recent demonstration that overexpressing 
OsAtlO (coincidentally also candidate 10 in our nomenclature) 
resulted in increased pCA almost certainly ester-linked to GAX 
(Bartley etal., 2013). This ester linkage is less common than FA 
in most tissues and is absent in wheat starchy endosperm. From 
all the evidence, we hypothesize that the activity responsible for 
GAX feruloylation will exist in some or all of the candidates 
1-5 and 8. 

The expression of candidate genes and amounts of bound FA 
and pCA were determined in tissues of developing B. distachyon 
plants as shown in Figure 2. 

The transcript abundances of the BAHD genes (denoted by the 
candidate numbers shown in Figure 1), were determined by quan- 
titative real-time RT-PCR (qRT-PCR) and the values are shown for 
the different genes, tissues, and developmental stages in Figure 3. 
The primer efficiencies were comparable for the different genes 
(Table 2), so differences in the abundances determined for the 
genes should reflect real differences in transcript abundances. On 
this basis, it can be concluded that the BAHD gene candidates 
which are most highly expressed in above-ground tissue samples 
from the vegetative phase are, in descending order of abundance, 
5, 9, 1, whereas in roots they are 1,9. In later stages of reproduc- 
tive development, they are 9, 7, 5, and 2p2, and in the spike at 
the advanced stage 7, 10, 5 (Figure 3). Therefore, BAHD genes 
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FIGURE 1 | Phylogenetic tree of clades within BAHD gene superfamily 
for rice and Brachypodium genes, which were identified as containing 
potential candidates for xylan feruloylation (Mitchell etal., 2007). 

Numbers are arbitrary candidate numbers based on rice genes; where two 
paralogous Brachypodium genes match one rice gene these are given the 



99 



B 



same number but with "p1 " or "p2" appended to the name. Genes within 
clade A are generally more highly expressed than those in clade B.The 
recently identified monolignol pCA transferase OsPMT (Withers etal., 2012) 
and OsAt10 (involved in pCA addition to xylan; Bartley etal., 2013) genes are 
indicated. 



from clade A are highly expressed in all tissues that were analyzed. 
Based on the phylogenetic relationship with OsPMT (Figure 1) 
and the assumption that BdBAHD 9 and 7 are responsible for the 
addition of coumaryl esters on lignin (i.e., PMT activity), and 
the others responsible for xylan feruloylation (with the excep- 
tion of 10 which is more likely responsible for coumaryl esters 
on xylan), it can be noted that at least one gene with each func- 
tion is highly expressed in every tissue that was analyzed, with the 
genes with PMT activity being relatively more highly expressed 
than other BAHD candidate genes in the reproductive growth 
phase. 

Determination of amounts of bound phenolic acids in plant 
samples is a measure of the amounts that are covalently linked to 
the cell wall fraction. Protocols using moderate alkaline treatment 
release only the ester-linked phenolics from cell walls, which are 
then separated and quantified by high-performance liquid chro- 
matography (HPLC). By far the most abundant bound phenolic 
acids are FA and pCA. FA is considered to be exclusively linked 
to GAX by ester bonds in grass cell walls, while the pCA is pre- 
dominantly ester-linked to lignin, with a much lower amount 
ester-linked to GAX. FA, unlike pCA, can oxidatively couple in 



planta to form dimers or higher-order oligomers. The amounts of 
bound FA monomer, major FA dimers and bound pCA are shown 
for different tissues in Figure 4. The amount of pCA reflects prin- 
cipally the amount of this phenolic acid which is ester-linked to 
lignin, and this fraction has been found to be highly correlated 
with the degree of lignification in grasses (Grabber etal., 2004). 
The values shown are therefore consistent with increased lignifica- 
tion during the reproductive phase, relative to GAX feruloylation, 
which occurs in both primary and secondary grass cell walls. The 
absolute amounts of FA and pCA in young shoots (early vegeta- 
tive phase, EVP) determined here are also comparable with those 
determined for seedlings at 8 days after germination (DAG) where 
FA and pCA were found to be ~400 and ~200 u,g g~ 1 , respectively 
(Christensen etal., 2010) and with a recent analysis of more tis- 
sues and stages (Rancour etal., 2012). A consistent feature of the 
data presented here (Figure 4) and by Rancour et al. (20 12) is that 
the pCA:FA ratio is higher in stems compared to leaves, which is 
consistent with the comparatively greater lignification in stems. 

Expression of BAHD candidates from both the PMT and non- 
PMT sub-clades was found in every tissue examined (Figure 3), 
as was substantial quantities of bound FA and pCA (Figure 4). 
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FIGURE 2 | Brachypodium developmental stages. (A) Early vegetative phase (EVP), bar = 1 cm; (B) late vegetative phase (LVP), bar = 
phase (TP), bar = 2 cm; (D) reproductive phase (RP), bar = 3 cm; (E) advanced phase (AP), bar = 2 cm. 



2 cm; (C) transition 



Simple correlations are not expected between the abundances of 
transcripts encoding enzymes in a biosynthetic pathway and the 
amount of end-product from that pathway. Transcript abundance 
measures the potential for enzyme synthesis only at a single time 
point whereas end-product accumulation occurs throughout a 
developmental process and may be limited by substrate availabil- 
ity as well as enzyme activity, which in turn depends on enzyme 
activation as well as enzyme amount. However, tissues which tend 
to have more bound FA compared to total bound pCA, might also 
be expected to have greater expression of BAHD candidates in the 
non-PMT sub-clade relative to those in the PMT sub-clade. We 
found this to be broadly the case, with shoots during vegetative 
phase and leaves in later phases, which have high bound total FA 
top pCA content (due to less lignification), also having greater 
ratios of non-PMT to PMT expression than internodes, roots, and 
spike (Figure 5). 

We have previously identified co-expression between BAHD 
candidates and xylan synthesis genes, compatible with a role in 
GAX feruloylation (Mitchell etal., 2007; Shewry etal, 2011), but 
given the new findings on these candidate genes (Withers etal, 
2012; Bartley etal., 2013), we decided to re-examine this. Public 
resources are not yet extensive for Brachypodium gene expression, 
but using the RiceFREND tool which employs a wide range of 



rice transcriptome experiments (Sato etal., 2013), we examined 
co-expression for the rice BAHD genes with rice genes putatively 
involved in xylan synthesis (Table 1). Nearly all the clade A genes 
had xylan synthetic genes in the top 1% of genes ranked by expres- 
sion correlation, but none of the clade B did. BAHD1, 5, and 
8 showed particularly close co-expression with genes involved 
in UDP-Xyl and UDP-Ara synthesis and GT61 genes, whereas 
OsPMT and similar genes (BAHD7 and BAHD9) had fewer closely 
co-expressed genes. If these genes in the PMT sub-clade are 
responsible for addition of pCA to lignin then some co-expression 
may be expected with any genes involved in secondary cell wall 
xylan synthesis. 

CONCLUSION 

The genes identified as candidates for being involved in feruloy- 
lation of GAX (Mitchell etal, 2007) are expressed in every tissue 
and at every developmental stage in the model grass Brachypodium 
as expected for a process which is required for every primary 
and secondary cell wall in the plant. The relative amounts of 
pCA and FA in these different tissues seem compatible with the 
expression of BAHD candidate genes, where these are divided 
into those putatively responsible for pCA linked to lignin (PMT) 
and FA linked to xylan (non-PMT). The evidence presented here 
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FIGURE 3 | Transcript abundance of candidate BAHD genes in developing 

Brachypodium tissues estimated by qRT-PCR. Error 

bars are ±SE from three biological replicates. One-way analysis 



of variance showed that the effect of tissue/stage on transcript 
abundance was significant at F probability < 0.001 for every gene except 
BAHD2p1, where it was significant at F probability = 0.002. 
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FIGURE 4 | Bound hydroxycinnamic acid content from developing 
Brachypodium tissues estimated by HPLC. Samples were early 
vegetative plants (EVP), late vegetative plants (LVP), transition phase 
leaves (TPL), transition phase internodes (TPI4, TPI3), transition phase 
roots (TPR), reproductive phase leaves (RPL), reproductive phase 
internodes (RPI4, RPI3, RPI2), advanced phase spike (APS), advanced 



phase leaves (APL), and advanced phase internodes (API4, API3, API2). 
FA dimers are the sum of the four major dimers described in Section 
"Materials and Methods." Error bars are ±SE from three biological 
replicates. One-way analysis of variance showed that the effect of 
tissue/stage on both total FA and pCA content was significant at F 
probability < 0.001. 
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FIGURE 5 | Ratios of expression of non-PMT BAHD to PMT BAHD and 
of bound FA to pCA. Expression ratio is the sum of values for genes 1-5 
and 8 (non-PMT) relative to the sum of values for genes 7 and 9 (PMT). FA 
to pCA ratio is the sum of FA monomer and dimers relative to pCA. Ratios 
shown are averaged over EVP LVP for shoots, TPL, RPL for leaves, TPI4, 
RPI4, RPI3, RPI2 for internodes; root is TPR and spike is APS. 



does not demonstrate function; for this more studies with tar- 
geted modification of the genes in grass species (as in Bartley 
etal., 2013) and in vitro characterization of pure proteins (as in 



Withers etal., 2012) are required. If their role in feruloylation is 
confirmed by such experiments, then as feruloylation is believed 
to be key to the exploitation of grass biomass for biofuel and ani- 
mal feed, these genes would represent important biotechnological 
targets. 

MATERIALS AND METHODS 
PLANT GROWTH 

Brachypodium plants Bd21 were grown in standard glass house 
conditions at 25°C. The whole aerial plant was harvested for the 
EVP (7 DAG) and late vegetative phase (LVP; 12 DAG). Plants 
were vernalized at 4°C for 2 weeks during the EVP; references to 
"DAG" do not include this period. In the transition (20 DAG) 
and reproductive (30 DAG) phases, the plants were separated into 
leaf and internode tissues, whereby internode 1 is the oldest and 
4 the upmost and youngest. In the transition phase, roots were 
also harvested. Spikelets, i.e., the developing seeds with the sur- 
rounding maternal tissues, were harvested in the advanced phase 
(50 DAG). Material was immediately frozen in liquid nitrogen and 
stored at -80°C. 

PHYL0GENETIC ANALYSIS 

Protein sequences for the whole BAHD superfamily were iden- 
tified as all rice, Arabidopsis, and Brachypodium sequences from 
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Table 2 | PCR primers. 



Oligo name 


Orientation 


Sequence (5'-3') 


Target 


Average primer efficiency 


prHMOl 


Sense 


GCCG C AC AAC AC C ATCATG 


qRT-PCR Amplicon Bd_Bahd_1 


1.74 


prHM02 


Antisense 


G G CTTTGATG AACTG CG CC 


qRT-PCR Amplicon Bd_Bahd_1 




prHM07 


Sense 


ACCTCATCCTCATGGCCCA 


qRT-PCR Amplicon Bd_Bahd_3p1 


1.745 


prHM08 


Antisense 


CGAAAACCAGGTGGCTGAAG 


qRT-PCR Amplicon Bd_Bahd_3p1 




prHM09 


Sense 


CCGGTAAGGTGGCCTCGT 


qRT-PCR Amplicon Bd_Bahd_4 


1.84 


prHMIO 


Antisense 


CGAACTCTGAAGGCAGCCG 


qRT-PCR Amplicon Bd_Bahd_4 




prHM11 


Sense 


CGTTCACCGCTTTCAACTTTG 


qRT-PCR Amplicon Bd_Bahd_5 


1.72 


prHM12 


Antisense 


TCGCAACCTGGTCTTTGACAC 


qRT-PCR Amplicon Bd_Bahd_5 




prHM19 


Sense 


CCGGTGCTAGCCCTGGAATA 


qRT-PCR Amplicon Bd_Bahd_10 


1.77 


prHM20 


Antisense 


TGCACGCGTTGTACTCCGA 


qRT-PCR Amplicon Bd_Bahd_10 




prHM33 


Sense 


CGCAAGACAATGACCGCTATG 


qRT-PCR Reference gene Bd_UBC18 


1.78 


prHM34 


Antisense 


CCAATCCGACGCCTCCTTATA 


qRT-PCR Reference gene Bd_UBC18 




prHM35 


Sense 


TGTTTGTGTCGGATTGGACGA 


qRT-PCR Amplicon Bd_Bahd_7 


1.69 


prHM36 


Antisense 


ACCGCCATGTAGTCCGCATAA 


qRT-PCR Amplicon Bd_Bahd_7 




prHM43 


Sense 


TTCTCGTATCACCCCTTCATGG 


qRT-PCR Amplicon Bd_Bahd_9 


1.75 


prHM44 


Antisense 


GGTGGTCTTCCTCCACACACAT 


qRT-PCR Amplicon Bd_Bahd_9 




prHM55 


Sense 


TGGCTTCTACGG C AACTG CTA 


qRT-PCR Amplicon Bd_Bahd_2p1 


1.79 


prHM56 


Antisense 


GCTTCCCGTCCTTGATGATCT 


qRT-PCR Amplicon Bd_Bahd_2p1 




prHM61 


Sense 


AAG C G G CTC G AGTAC CTG 


qRT-PCR Amplicon Bd_Bahd_2p2 


1.82 


prHM62 


Antisense 


GCCATTGTTGCTGGAGTTGT 


qRT-PCR Amplicon Bd_Bahd_2p2 




pin mod 


Sense 


PP ACCTCTCCTTCCCC ATP 


cjn i-rLn Ampncon Da_Dana_o 


1 7Ct 

I. /a 


prHM66 


Antisense 


CGCATGATGTAGTAGCAGTTGCC 


qRT-PCR Amplicon Bd_Bahd_8 




prHM69 


Sense 


ATCCCGCCATCCAACATCTAC 


qRT-PCR Amplicon Bd_Bahd_12 


1.79 


prHM70 


Antisense 


CG GATCTG G CCTTG ATGTTGT 


qRT-PCR Amplicon Bd_Bahd_12 




prHM75 


Sense 


TTCAACAGCATGGATGGCC 


qRT-PCR Reference gene Bd_SDH 


1.76 


prHM76 


Antisense 


ATCTTCGGTTGCAGAGCTCCT 


qRT-PCR Reference gene Bd_SDH 





Phytozome 1 containing the PFAM domain PF02458. An initial 
tree was generated and all sequences contained within the candi- 
date group identified in Mitchell et al. (2007) were aligned with 
the MUSCLE algorithm (Edgar, 2004). Gapped columns were 
removed, followed by phylogeny analysis of aligned sequences in 
the PhyML package (Guindon etal, 2005) using the Whelan and 
Goldman (2001) model. An initial run optimized the gamma and 
invariant proportion parameters; these were then held constant for 
100 runs for bootstrap non-parametric analysis and the maximum 
likelihood tree is presented. 

DETERMINATION OF TRANSCRIPT ABUNDANCE 

RNA was extracted using a cetyltrimethylammonium bromide 
(CTAB) method following Chang etal. (1993). Quantitative 
PCR was performed as in Pellny etal. (2008) on an Applied 
Biosystems 7500 real-time PCR system 2 using SYBR green Jump- 
Start Kit (Sigma-Aldrich 3 ) following the manufacturer's instruc- 
tions. The expression of the genes of interest was normalized 



1 www.phytozome.org 

2 www.appliedbiosystems.com 

3 www.sigmaaldrich.com 



with two endogenous controls ubiquitin-conjugating enzyme 
18 (UBC18) and succinate dehydrogenase (SDH). Putative 
house-keeping genes glyceraldehyde-3-phosphate dehydrogenase 
(GAPDH), S-adenosylmethionine decarboxylase (SAMDC), and 
peptide deformylase 2 (PDF2) were also tested but found less stable 
over the different tissue and time points. Expression calculations 
were performed as suggested by Rieu and Powers (2009). Indi- 
vidual amplification efficiencies were established with LinRegPCR 
using a window-of-linearity and the relative quantities (RQs) for 
the target genes were calculated. This was normalized with the 
geometrical mean of the RQs of the control genes. Primers used 
and average efficiency of pairs are listed in Table 2. 

DETERMINATION OF BOUND PHENOLIC CONTENT 

Cell wall bound phenolics were released by alkaline hydrolysis of 
alcohol insoluble residues ( AIRs) from samples (20 mg) of freeze- 
dried, ground tissue and extraction into ethyl acetate as previously 
described (Pellny etal, 2012). 3,5-Dichloro-4-hydroxybenzoic 
acid (20 |iL at 1.5 mg mL _1 ) was added to AIR of all samples, 
prior to alkaline hydrolysis, as an internal standard. Samples were 
dissolved in 1 mL 50% methanol:2% acetic acid (v:v) and bound 
phenolic acids from 40 u.L of extract separated by HPLC on a 
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Shimadzu Prominence high-performance liquid chromatograph 
as described by Waldron etal. (1996) but using a binary gradient 
pump system. Quantitation of FA and pCA was by integration 
of peak areas at 280 nm with reference to calibrations made using 
known amounts of pure compounds. Peaks of the major FA dimers 
were identified by comparison of retention times with pure stan- 
dards kindly supplied by Professor John Ralph (Lu et al., 2012; 5-5, 
8-0-4, 8-5-benzofuran) or by comparison of spectrum with that 
in Waldron etal. (1996) (8-5). Dimer quantitation was achieved 
relative to FA monomer using areas of these peaks and the response 
factors for dimers and FA monomer in Waldron etal. (1996). All 
samples were extracted and analyzed in triplicate. 
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